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One important question in relativistic heavy ion collisions is if hadrons, specifically anti-hyperons, 
are in equilibrium before thermal freezeout because strangeness enhancement has long been pointed 
to as a signature for Quark Gluon Plasma. Because anti-baryons have long equilibration times in the 
hadron gas phase it has been suggested that they are "born" into equilibrium. However, Hagedorn 
states, massive resonances, which are thought to appear near the critical temperature, contribute 
to fast chemical equilibration times for a hadron gas by providing extra degrees of freedom. Here 
£ — , we use master equations to describe the interplay between Hagedorn resonances, pions, and baryon 

anti-baryon pairs as they equilibrate over time and observe if the baryons and anti-baryons are fully 
equilibrated within the fireball. 
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I. INTRODUCTION 

m . 

It has been suggested that one of the main signatures for the Quark Gluon Plasma (QGP), which is a state of 
matter that consists of deconfined quarks and gluons, could be the equilibration of anti-hyperons and multi-strange 
baryons [![. The idea was that within a hadron gas there would not be enough time for these particles to reach 
chemical equilibrium and thus they should be "born" into equilibrium after the QGP phase transition 0,0]. Indeed, 
\ experimentally we know that the particle abundancies reach chemical equilibration close to the phase transition 4]. 
Multi-mesonic reactions within the standard hadron gas model, which can explain abundancies at SPS energies @, 6], 
cannot account for the short chemical equilibration times and high abundancies of anti-baryons at RHIC [7|, [8| (but 
also see @). In this paper we focus on Hagedorn resonances (heavy resonances that have an exponential mass 
spectrum that appear near the critical temperature) that drive multi-pionic reactions and also produce baryon anti- 
baryon pairs as was suggested in Ref. Once we include the Hagedorn resonances we see that it is possible for 
the anti-baryons (future work will discuss anti-hyperons and multi-strange baryons) to quickly equilibrate within the 
T ! fireball just "below" the phase transition. 

Originally, (anti-)strangeness enhancement at CERN-SPS energies in comparison to pp-data, which was primarily 
observed in anti-hyperons and multi-strange baryons, was thought of as a signature for QGP. Using binary strangeness 
■ production reactions such as 



> 
OS 



o 



X 



or binary strangeness exchange reactions 



p^K + K (1) 



K + p <-> 7r + A (2) 



it was concluded that it took far too long for chemical equilibrium to be reached within the hadron gas phase 
Thus, it was proposed that the QGP was already observed at SPS because strange quarks can be produced more 
abundantly by gluon fusion, which would then account for strangeness enhancement in hadrons following hadronization 
and rescattering of strange quarks [l[ . 

On the other hand, to explain secondary production of anti-hyperons the idea was suggested that strangeness 
enhancement could be explained using multi-mesonic reactions such as 

p + N <-> nn (3) 
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for the anti-protons § and for the anti-hyperons 



£,A + A <-> mr + AT 
S + A «-> mr + 2K 
tt + N <-> mr + 3AT 



found in Ref. Q . The anti-hyperons can be rewritten into the general equation 

Y + N «-> n-n- + nyif. 



(4) 



(5) 



The arrows in Eqs. ([3]), (J4j) signify the equal probability that a decay can occur in either direction otherwise known 
as detailed balance. 

The time scale of a standard hadron gas at SPS can be then estimated using the multi mesonic reactions in Eq. 
©, Q. Using 



c jvy ~ &Np 



50 mb 



the chemical equilibration time is then 



7y = (T-5 



((cr NY v NY ))pi 



1 



fm 



(6) 



(7) 



where w po to 2po, which is typical for SPS 0, HI- Therefore, the time scale given in Eq. ([7]) is short enough to 
account for chemical equilibration within the cooling hadronic fireball at SPS. 

If we apply our same understanding of the hadron gas phase to RHIC temperatures our time scales are much longer. 



The equilibrium rate of at RHIC at T = 170 MeV is Y c £ 



the baryon density is N'g — « 0.04 fm 3 , which leads to a time scale of tq = 



( a nB v ns)^B where the cross section is a ~ 30 mb and 



10 — . However, considering 



that the fireball's time scale in the hadronic stage is r < 4 — , a standard hadron gas could not explain the apparent 
chemical equilibration observed in baryons within the fireball. Moreover, these results were also backed up in Ref. Q 
where using a fluctuation-dissipation theorem it was found that the equilibration time of baryons and anti-baryons, 
r « 10 [— 1 , at RHIC temperatures. In the 5% most central Au-Au collisions the baryon anti-baryon production is 
roughly three times lower than the measured experimental values if it starts out of equilibrium (specifically at zero 
for reactions of type (O and (H])) shown in Fig. |T]), taken from Ref. [8J, for A production. Because of the apparent 
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FIG. 1: Anti-lambda production in the most central collisions at RHIC. The bottom three dashed lines start with no anti- 
lambdas at T — 180 MeV whereas the solid lines assume that the anti-lambdas begin in equilibrium. The shaded area shows 
the experimental results taken from Ref. |8j. The set of lines shows the difference for varying coupling coefficients. 

differences in the equilibration times some have suggested that the hadrons are "born" into equilibrium i.e. the system 
is already in a chemically frozen out state at the end of the phase transition [3J] . 

In this paper we use Hagedorn states to provide the extra degrees of freedom needed to match experimental results 
when the anti-baryons begin out of equilibrium. Baryon anti-baryon production develops according to the possible 
reaction 



7i7T <-> HS <-> nn + BB, 



(8) 
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which is initially discussed in Ref . [10( . The Hagedorn states arethought as massive resonances with short time scales 
and only contribute near the critical temperature. They can then catalyze rapid equilibration of baryons and anti- 
baryons near T c . The results obtained here lead us to believe that the baryon anti-baryon pairs have sufficient time 
to equilibrate within the fireball close to the phase boundary. 



II. HAGEDORN STATES 

In the 1960's Hagedorn found a fit for an experimentally growing mass spectrum [TTj . p, which is described as 

pin m 
F(m)e T H dm, 

F{m) = -. (9) 

(m 2 + m 2 ,) 4 

where the minimum mass is mo — 500 MeV, the maximum mass is M = 7 GeV, the Hagedorn temperature is 
T H = 180 MeV, which fits within Lattice QCD predictions [13], and A = 0.5 is a free parameter. Because of the 
exponential growth, near the Hagedorn temperature the Hagedorn states can account for the extra degrees of freedom 
needed to match experimental values. A is then chosen by looking at the energy density and trying take into account 
the extra degrees of freedom needed. Here we are considering only mesonic, non-strange Hagedorn states with masses 
between 2 GeV and 7 GeV. 

To describe the dynamical behaviour of the Hagedorn states we use rate equations. Rate equations have both loss 
and gain terms and have the basic form 

dn 

— = — loss + gain. (10) 
dt 

We assume a system with no net baryon density i.e. N B = N B = N BB . This should approximately be the case at 
RHIC at midrapidity. For the reaction HS <-> nil + BB the behaviour of the density of the i th Hagedorn resonance 
Nru\, pions iV-Tr, and baryon anti-baryon pairs N BB is described by the following set of equation: 



dt 



n 

rr BB ^if B (T)(N^Ni 



i,BB'-"-i,(n) V- 1 JV-'nJ BB 

dN^ 



^^TftnB^ (N m - 5R(T)(JV ff ) 



^BB 



dN BB 
dt 



i 

-Y, T Vbb (N 2 BB N^M h{n) (T) - N m ) (11) 



where T is the decay width, _B.^ n7r represents the branching ratios, and (n) is the average number of pions that 
each Hagedorn state decays into when a baryon anti-baryon pair is included. We also have two separate detailed 

balance factors, !ft(T) = for the decay HS <-> nir and ^|™^ B ^(T) = jy^nvW 5 f° r tne decay HS <-> nir + BB. 

The detailed balance factors ensure that detailed balance is maintained in equilibrium. They are also temperature 
dependent because the equilibrium values of the density are dependent on the temperature. 
The branching ratios are described by a gaussian distribution 

1 W) 2 , , 

Si-™**— r=e (12) 

<7V27T 



where (n) = 0.6 + 0.3^- 5— 16 is the average pion number that each Hagedorn state decays into and a = 0.26^^- 

is the width of the distribution [l3j. The decay width = 0.17toj — 88, which has the range Tj = 250 MeV to 1100 
MeV, is a linear fit extrapolated from the data from the particle data group [T3|. In the future we will also use a 
microcanonical model to find the branching ratios as shown in Ref. [l5j. 
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III. RESULTS 



We have solved the rate equations in Eq. (fTTj) considering several different initial conditions for a statistical system. 
At first we take the simplest example and observe only the decay HS <-> nir when the pions are held in equilibrium. 
We then consider a resonance bath and allow the pions to equilibrate. We also allow both the pions and the resonances 
to develop until they reach equilibrium. Then we include baryon anti-baryon pairs into our decay HS — > wk + BB 
where at first the pions are held in equilibrium, then the Hagedorn states, and also both the Hagedorn resonances and 
the pions are held in equilibrium while the baryon anti-baryon pairs are allowed to equilibrate. Finally, we consider 
the case when all the constituents are allowed to equilibrate simultaneously. 



A. Case 1 HS <-> mr: Pions held in Equilibrium 

Initially, when the Hagedorn states decay only into multiple pions, HS <-> nn we assume that the pions start 
in equilibrium and the Hagedorn resonances start at zero, which could be an approximation for a physical system 
immediately following hadronization when the correlation lengths are very short. Then we can make an estimate 
for the time scale with the inverse of the decay width T = -. The equilibration time estimate is then between 
t w 0.2 — and 0.5 — . In Fig. @ the result of the rate equations is shown, which matches the time estimates well. 




FIG. 2: Hagedorn resonances for various masses at T — 170 MeV. 



From Fig. ([2]) we see that the Hagedorn states equilibrate quickly (in comparison to typical expansion times). Hence, 
they should be at chemical equilibrium as long as the pions are at chemical equilibrium. 



B. Case 1 HS <-> nn: Hagedorn Resonances held in Equilibrium 

On the other hand, if the Hagedorn resonances are held in equilibrium and treated like a resonance bath then we 
can define a new effective production rate for pions: 

N eq 

i 

which gives the range for the equilibration time r « 0.02 — to 0.25 — . Here the pions start at zero, which could be 
an approximation for a physical system immediately following hadronization when the correlation lengths are very 
long. Again the time scale can be compared with the results from the rate equations, which is shown in Fig. ([3]). 
Indeed, the pion become populated very quickly towards chemical equilibrium for all temperatures. 
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Case 2: HS in Eq. 

MM=7 GeV T<,=180 MeV 




— T=160 MeV T=0.269 fm/c 
■■ T=165 MeV T=0.138 fm/c 

- ■ T=170 MeV T=0.069 fm/c 
■- T=175 MeV T=0.034 fm/c 

T=180 MeV T=0.017 fm/c 



0.5 

Time (fm/c) 



0.75 



FIG. 3: Number density of pions as they go to equilibrium while the Hagedorn resonances are held in equilibrium. 

C. Case 3 HS «-» nn: Hagedorn Resonances and Pions are both out of Equilibrium 

The most interesting case is when both the Hagedorn states and the pions are out of equilibrium. The coupled 
network of rate equations is 



dN 



R{i) 



(It 



At 



oo oo 

= E E n ot nB^ nn (N m - »(T)(JV T )»). 



(14) 



i n—2 



The right-hand side of the rate equations goes to zero when the particles are in equilibrium. In Eq. (|14[) the resonances 
and pions reach a quasi-equilibrium configuration before full equilibrium is reached. In the quasi-equilibrium state 

the right-hand side nears zero and thus slows down the equilibration time, thus, N R ^ w ^R{i) 



(n) 



and Ntt 



N% q [ ) . In Fig. ((4]) quasi-equilibrium is quickly reached on the time scales from case 1 and case 2. 
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Case 3: 7l's: HS: 2*Eq 
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Case 2: tc's: HS: 2*Eq 



eVT K =18 
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FIG. 4: Here the resonances are at twice the equilibrium values and the pions start at zero. 
The quasi-equilibrium can be understood by looking at the number of effective pions in the system i.e. 

N* = N« + '%2(n)N R (15) 

R 

where (n) is the average number of pions that the Hagedorn resonance Nn will decay into. Essentially in this particular 
example, there are too many effectives pions in the system and they must be killed off, which can account for the 
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longer time scales. The direct pions are practically already in equilibrium after t > 0.5 — , especially at higher 
temperatures, however, the Hagcdorn states are clearly overpopulated, which implies that the effective number of 
pions are also overpopulated. Thus, the long time scale consists of reactions such as such as nir — > HS — > mir 
where n > m, which kills off the effective pions in the system. Moreover, the heavier Hagedorn resonances are more 
overpopulated because they have a higher probability to decay into more pions than lighter resonances. A more in 
depth look into the time scales of the effective pion number will be discussed in an upcoming paper. 



D. HS <-> nir + BB 

Finally, we want to determine the production of baryon anti-baryon pairs close to the phase transition. The 
baryon anti-baryon pairs are produced through the reaction found in Eq. (jHJ) and these particular branching ratios 
are determined using a microcanonical model. The average pion number for each corresponding mass is calculated by 
making a fit to the multiplicities in the microcanonical model found in Ref. [l5l ]. 




FIG. 5: Various baryon multiplicities, which are used to determine the total baryon density in our system, for the decay 
HS^nw + BB^,^. 

For baryon anti-baryon pairs we need to also reconfigure the decay widths, which is also done with the baryon 
multiplicities from Ref. ElJ. We assume that for every baryon there is a corresponding anti-baryon and then we 
add up the multiplicities for the three dominat baryons: the proton, the neutron and lambda, to determine the total 
baryon multiplicity (B) as shown in Fig. ([5]) We can then estimate the relative decay width for the baryon anti-baryon 
decay, r* ^ = {B)T\°*. The total baryon multiplicity varies between 0.2 and 0.4, which then gives a decay width 

between Tf^ = 50 MeV and 450 MeV [nj. 

We start with case 1 when both the Hagedorn resonances and the pions can be held in equilibrium, while the baryon 
anti-baryon pairs are driven to equilibrium as shown in case 1 in Fig. ([6]). Case 2 is again when the pions are held in 



Case 1: TC's and HS in Eq. 



MM=7 („'Y ■] =!-!! Mu\ 




BB T= 1 60 Me V T= 1 .062 fm/c 
BB T-165 MeV 1=0.696 fm/c 
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Case 2: 7C's in Eq. 

MM=7 GeV T =180 MeV 
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Case 3: HS in Eq. 

MM=7GcV T =180 MeV 
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BBT=180Mev 
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FIG. 6: Case 1: Solution of the baryon anti-baryon rate equation when the pions and Hagedorn states are held in equilibrium. 
Case 2: Baryon anti-baryon pairs when the pions are held in equilibrium. Case 3: Baryon anti-baryon pairs when the Hagedorn 
states are held in equilibrium. 
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equilibrium (and the resonances begin at zero) and, case 3 is when the Hagedorn resonances are held in equilibrium 
(and the pions are started at zero). As before we can make a time scale estimate, this time for the baryon anti-baryon 
pairs. Our effective chemical equilibration for the baryon anti-baryon pairs is then 

/ N eq \ 

T eff _ T tot _rm (16) 

i,BB 2—i i,BB I N eq_ I ' V lu ) 

i \ BB / 

which gives equilibration times r = 0.2 — to 1 — . The results of the pions being held in equilibrium are shown in 
case 2 in Fig. ((6]) and the results for the Hagedorn states held in equilibrium are shown in case 3 in Fig. ([6]). 

For case 2 the baryon anti-baryon pairs take slightly longer to equilibrate than their time scale estimate because 
quasi-equilibrium is reached with the Hagedorn states, which occurs around r ~ 1 ^p. Case 3 is not as affected by a 
quasi-equilibrium state because the pions equilibrate so quickly that they do not affect the equilibration time of the 
baryon anti-baryon pairs. Comparing the graph of baryon anti-baryon pairs in case 1 to that in case 2 in Fig. ([6]) 
we clearly see that they are still almost identical. The reason is that the baryon anti-baryon pairs are not affected 
by the pions because the pions equilibrate almost immediately and thus the approximation that the pions are held in 
equilibrium can be made. 

When everything is allowed to develop out of equilibrium, quasi-equilibrium is reached on a time scale of the 
previous estimated equilibration times. The results are given in Fig. J7]| Although we only see a small deviation in 
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Case 2: HS=2*Eq jt=0 BB=0 
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Case 4: HS=2*Eq 71=0 BB=0 
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FIG. 7: Solution of the rate equations when the resonances, pions, and baryon anti-baryon pairs are allowed to change over 
time. 

the pions from equilibrium in Fig. (|7|) that minor deviation drastically affects the resonances because even the lightest 
resonance, whose mass is M = 2 GeV, decays on average into (n) 5 pions. In Fig. ([7]) we clearly see that quasi- 
equilibrium is reached quickly i.e. T quas i- eq < 1 — . After quasi-equilibrium is reached the remaining constituents 
(especially the resonances) slowly reach equilibrium. What we can get from Fig. ([7]) is that the pions and the baryon 
anti-baryon pairs quickly equilibrate, especially for higher temperatures. In the graph of the baryon anti-baron pairs 
we see that at T — 1 80 MeV the baryon anti-baryon pairs quickly populate near chemical equilibrium and while they 
are not in complete equilibrium are quite close to it. Hence, for temperatures between T = 180 MeV and 160 MeV, 
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the baryon anti-baryon pairs are populated with equilibration times faster or equal to 2 isi. This constitutes a very 
promising finding. 



IV. CONCLUSIONS 



Our preliminary results and time scale estimates indicate that baryon anti-baryon pairs can be "born" out of 
equilibrium after hadronization and then equilibrated by the subsequent population and decay of Hagedorn states. 
Cases 1-3 for the decay HS «-> nir + BB clearly show quick equilibration times for baryon anti-baryon pairs between 
temperatures of 180 MeV and 160 MeV i.e. slightly below the phase transition. When all the particles started out of 
equilibration the baryon anti-baryon pairs quickly neared equilibrium although a minor deviation was still observed 
from chemical equilibration. Afterwards, due to affects from the need to kill of the number of effective pions in the 
system, longer time scales were observed when the pions, Hagedorn states and baryon anti-baryon pairs were out of 
equilibrium. However, the pions and baryon anti-baryon pairs were quickly populated near T c and remained close to 
their equilibrium values even when it took longer for the Hagedorn states to reach chemical equilibrium. Since the 
Hagedorn states only contribute near T c these appear to be acceptable results. 

In an upcoming publication we will delve more thoroughly into the effects of the quasi-equilibrated state seen in both 
Fig. (|3]), and we also will discuss its effects on the fireball as it cools over time due to a Bjorken expansion. Thus, 
we expect to see the baryon anti-baryon pairs quickly reach equilibrium and then they will not be produced further 
when the system is cooled. Moreover, we will vary our choices of initial conditions and discuss their implications. 
Finally, a non-zero strangeness will be included specifically in the baryon anti-baryon pairs so that we can study the 
equilibration times of anti-hyperons and multi-strange baryons. The continuation of our work is promising in order 
to explain results for strangeness in baryons and anti-baryons found at RHIC. 
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